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Abstract 

The path integral of the relativistic Coulomb system is solved, and the wave 
functions are extracted from the resulting amplitude. 
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1) While the path integral of the nonrelativistic Coulomb system has been solved some 15 



years ago [HI and further discussed by many authors B - so that it has become textbook 



material ||15|| , the relativistic problem has remained open. The purpose of this note is to fill 
this gap. 

2) Consider first a relativistic free particle of mass M. If x(X) describes its orbit in D 
spacetime dimensions in terms of some parameter A, the classical action reads 



Mc 



(1) 



where c is the light velocity. This action cannot be used to set up a path integral for the 
time evolution amplitude since it would not yield the well-known Green function of a Klein- 
Gordon field. An action which serves this purpose can be constructed with the help of an 



auxiliary fluctuating variable p(A) and reads [16| 



A 



dX 



M 



Mc 2 



P(A) 



_2p(A) vv ' 2 

Classically, this action coincides with the original A c \ since it is extremal for 

p(A) = y/^Xj/c. 

Inserting this back into (§) we see that A reduces to A c \. 

The action ^4 c i is invariant under arbitrary reparametrizations 

A -/(A). 

The action A shares this invariance, if p(A) is simultaneously transformed as 

P - p/f- 



(2) 



(3) 



(4) 



(5) 



The action A has the advantage of being quadratic in the orbital variable x(A). If 
the physical time is analytically continued to imaginary values so that the metric becomes 
euclidean, the action looks like that of a nonrelativistic particle moving as a function of a 
pseudotime A through a D-dimensional euclidean spacetime, with a mass depending on A. 



3) To set up a path integral, the action has to be pseudotime-sliced, say at Ao = 
A a , Ai, . . . , Atv+i = A&. If e n = A ra — A n _i denotes the thickness of the nth slice, the sliced 
action reads 



N+l 
n=l 



M . N2 Mc 2 
[Ax n ) + — — e n p n 



(6) 



where p n = p(A ra ), Ax n = x n — x n _i, and x n = x(X n ). A path integral f(V D x/y/p D )e A ^ h 
may be defined as the limit iV —>■ oo of the product of integrals 



-D 



N 

n 



2nhe b p b /M n= 



d D 



-D 



2irhe n p n /M 



(7) 



This can immediately be evaluated, yielding 



-D 



2M/Mc 



cxp 



Mc (x b - x a ) 2 Mc 



(8) 



where the quantity 



has the continuum limit 



N+l 

L = C^ tnPn 
n=l 



L = c d\p(X). 



(9) 



(10) 



Classically, this is the reparametrization invariant length of a path, as is obvious after 
inserting (JJ). 

If the amplitude (§) is multiplied by Ac/2, where Ac = h/Mc is the Compton wavelength 
of the particle, an integral over all positive L yields the correct Klein-Gordon amplitude 



1 / Mc 



D/2-l 



Kd/2-! (McVx^/h) 



(11) 



(2tt)^/ 2 {hV^J 
where K v (z) is the modified Bessel function. 

The result does not depend on the choice of p(A), this being a manifestation of the 
reparametrization invariance. We may therefore write the continuum version of the path 
integral for the relativistic free particle as 



I r°° t r T) D r 



-A 



(12) 



2 -Jo J ' L ' J J J? 
where $[p] denotes a convenient gauge-fixing functional, for instance $[p] = S[p — 1] which 
fixes p(A) to unity everywhere. 

To understand the factor l/^fp D in the measure of (|12|), we make use of the canonical 
form of the action (0), 



A[p,x] = I dX 

J A, 



After pseudotime slicing, it reads 



. , p(A)p 2 Mc 2 . 



(13) 



JV+l 
n=l L 



+ 



2M 2 

At a fixed p(A), the path integral has then the usual canonical measure: 



(14) 



N 

n 

n=l 



-, N+l 

n 



n=l L 



(2nh) 



D 



-A N [p,x]/h 



{2kK) d 

By integrating out the momenta, we obtain (0) with the action (|6|). 

4) The fixed-energy amplitude is related to (O) by a Laplace transformation: 



(15) 



poo 

(x 6 |x a ) E = -i / dx b e E ^-*°V n (x b \x a ) 



(16) 



where x° denotes the temporal component and x the purely spatial part of the D-dimensional 
vector x. The poles and cut of (xb|x a ) E along the energy axis contain all information on the 
bound and continuous eigenstates of the system. The fixed-energy amplitude has the path 
integral representation 

V D - l x 



J y/p 



-Ae/K 



(17) 



with the action 



A E = dX 



M /n E 2 Mc 2 ... 

-x (A) - p(A)— — + -r-p(X) 



(18) 



2p(A) v ' rv J 2Mc 2 ' 2 
This is seen by writing the temporal part of the sliced D-dimensional action (|6[) in the 
canonical form fllil) . By integrating out the temporal coordinates x^ in (|T5|) , we obtain 



iV (^-functions. These remove N integrals over the momentum variables leaving only a 
single integral over a common p°. The Laplace transform (0), finally, eliminates also this 
integral, making p 4 equal to —iE/c. In the continuum limit, we thus obtain the action (|T8|). 
5) The path integral ( |TTD forms the basis for studying relativistic particles in an external 
time-independent potential V(x). This is introduced into the path integral ( |T7| ) by simply 
substituting the energy E by E — V(x). 

For an attractive Coulomb potential in D — 1 = 3 spatial dimensions, the above substi- 
tution changes the second term in the action ( |18|) to 

The associated path integral is calculated with the help of a Duru-Kleinert transformation 
IJ as follows: 

First, we increase the three-dimensional configuration space in a trivial way by a dummy 
forth component x 4 (as in the nonrelativistic case). The additional variable x A is eliminated 
at the end by an integral / <ir 4 /|x a | = / dj a (see Eqs. (13.114) and (13.121) in Ref. [0). 
Then we perform a Kustaanheimo-Stiefel transformation dx^ = 2A(u) fl u du l/ (Eq. (13.100) in 
Ref. yi5| ). This changes x'^ 2 into 4u 2 u' 2 , with the vector symbol indicating the four-vector 
nature of u. The transformed action reads 



r*t> r4Mu 2 ^ /2 . x . p(A) 



P 4 

n 2 



(M 2 c 4 -^)r-2Ee 2 



(20) 



2p(A) v ' 2Mc 2 u 2 

We now choose the gauge p(A) = 1, and go from A to a new parameter s via the time 
transformation dX = fds with / = u 2 . This leads to the Duru-Kleinert-transformed action 

s b r4Af , 0/ x l 



A 



DK 

E 



It describes a particle of mass fi = 4M moving as a function of the "pseudotime" s in a 
harmonic oscillator potential of frequency 

lj = ^VM2 C 4 - E \ (22) 
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The oscillator possesses an additional attractive potential — e 4 /2Mc 2 M 2 which is conveniently 
parametrized in the form of a centrifugal barrier 

Kxtra = ft 2 |Nt, (23) 

whose squared angular momentum has the negative value 

4 tra = -4a 2 , (24) 

where a denotes the feinstructure constant a = e 2 /he ~ 1/137. In addition, there is also a 
trivial constant potential 

Konst = -Jj^ e ■ ( 25 ) 

If we ignore, for a moment, the centrifugal barrier V extra , the solution of the path integral 



can immediately be written down (compare Eq. (13.121) in Ref. [ |T5|j ) 

(x b \x a ) E = - t —-j o dLe eEL / Mch J Q d la (u b L\u a 0), (26) 

where (v,bL\u a 0) is the time evolution amplitude of the four-dimensional harmonic oscillator. 

There are no time-slicing corrections for the same reason as in the three-dimensional 
case. This is ensured by the affine connection of the Kustaanheimo-Stiefel transformation 
satisfying 

T/ A = <Te 4 W = (27) 



[see the discussion in Section 13.6 of Ref. |L5 
A 7 a -integration leads to 

h Mk f 1 , Q~ v _ / r — N 



-K- (r b + r Q 

1 - Q 



(28) 



x exp 
with the variable 

p = e~ 2ujL , (29) 



and the parameters 

e 2 E a 



v 



2uhMc 2 Jm 2 c 4 /E 2 - 1 



2h he hcu K J 



We now use the well-known expansion 



2 00 

I o (zcos(0/2)) = ~£(2Z + l)fl(cos0)I a+ i(s) (31) 

2 1=0 



and obtain the partial wave decomposition 



1 °° 21 + 1 

(x 6 |x a )£ = y2{r b \r a ) E j— Pi{cos8) 

r b r a ^ 4tt 

1 oo I 

= E( r b\r a ) E ,i £ *W*&)*z™(*«)> ( 32 ) 

with the usual notation for Legendre polynomials and spherical harmonics. The radial 
amplitude is, therefore, 

% 2M f°° 1 

(ThWaiE i — — i \/rhr a —— / dy — : — e 2uy (33) 

x exp [-Kcothy(r 6 + r a )\ I 2 i+i f 2k V r & r a ^ ■ 

At this place, the additional centrifugal barrier (^) is incorporated via the replacement 



21 + 1 - 21 + 1 = J(2Z + l) 2 + 4tra- (34) 



(as in Eqs. (8.146) and (14.237) in Ref. ||15|| ). The integration over y yields 

(r b \r a ) El = — »— ^— r ^~1 + 1 + ^ Wr,^ (2rcr b ) M 1/2 (2«r a ) (35) 
1 1 ; ' 2Mck (2Z + 1)! " ,i+1/2 ^+1/2 v a) y , 



(compare Eq. (9.52) in Ref. |5|). 

This fixed-energy amplitude has poles in the Gamma function whenever v — I — 1 = 
0,1,2,... . They determine the bound-state energies of the Coulomb system. Subsequent 
formulas can be simplified by introducing the small positive /-dependent parameter 

5l = l_l = l + 1/2 - J (i + 1/2)2 _ a 2 w + (36 ) 



With it, the pole positions are given by v = hi = n — Si, with n = I + 1, 1 + 2, 1 + 3, . . . , and 
the bound state energies become: 



E 



id 



±Mc 2 
±Mc 2 



1 + 



-1/2 



(n - 5 l )\ 



1 3 



(37) 



Note the appearance of the plus-minus sign as a characteristic property of energies in rel- 
ativistic quantum mechanics. A correct interpretation of the negative energies as positive 
energies of antiparticles is straightforward within quantum field theory; it will not be dis- 
cussed here. 

To find the wave functions, we approximate near the poles v ~ hi: 



r(-i/ + z + i; 
i 



n T \ v — h\ 



2 K 2 k 2 ( E \ z 2Mc 2 



u-ni 

K PS 



hi 2M 
E 1 1 



\Mc 2 ) 



E 2 -E^ 



(38) 



Mc 2 an hi ' 

with the radial quantum number n r = n — I — 1 . In analogy with a corresponding nonrela 



tivistic equation (Eq. (13.203) in Ref. fl5fl), the last equation can be rewritten as 

1 1 



K 



a H v 



(39) 



where 



Mc 2 



Off = CLh~ 



E 



(40) 



denotes a modified energy-dependent Bohr radius. Instead of being 1/a w 137 times the 
Compton wave length of the electron h/Mc, the modified Bohr radius which sets the length 
scale of relativistic bound states involves the energy E instead of the rest energy Mc 2 . 

With the above parameters, the positive-energy poles in the Gamma function can be 
written as 



-iT(-v + l + 1 



-) n - 1 / E \ 2 2Mc 2 ih 



M 



Tik nfn r \ a H \Mc 2 J E 2 - E 2 nl 



(41) 



Using this behavior and a property of the Whittaker functions (see Eq. (9.73) in Ref. ||15||), we 
write the contribution of the bound states to the spectral representation of the fixed-energy 
amplitude as 



(42) 



n=l+l X±VJ -^ ' ^ ^nl 

A comparison between the pole terms in ([35]) and fl42 ) renders the radial wave functions 



Rniir) 



affn (2l + l)\\(n-l-l)\ 



x(2r/nia H ) 



l+l -r/ni&H 



M(-n + 1 + 1,21 + 2, 2r/hia H ) 



(43) 



(n-l-l)\ _ r/ 



(2r/n l ~aH) l+1 L 2 ^l 1 (2r/n l ~aH) 



~"¥ni\ (n + iy 



The properly normalized total wave functions are 



1 



^nZm(x) = -R n i(r)Yi m (x). 
r 



(44) 



The continuous wave functions are obtained in the same way as in the non-relativistic 
case (see Eqs. (13.211)-(13.219) in Ref. [fTJ). 

This concludes the solution of the path integral of the relativistic Coulomb system. 
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